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Abstract 


It is well known that there are several processes to manufacture composite materials, 
a large part of which consist in the infiltration of a liquid (matrix) through a porous 
medium (reinforcement). To perform these processes, both thermodynamics (wet- 
ting) and kinetics (Navier-Stokes) must be considered if a good quality composite 
material is sought. Although wetting and the laws that govern it have been well known 
for over 200 years, dating back to the original works of Young and Laplace, this is not 
the case with the Navier-Stokes equation, which remains so far unsolved. Although 
the Navier-Stokes equation, which describes the motion of a fluid, has been solved for 
many particular cases, such as the motion of a fluid through a pipe, which has resulted 
in the well-known Poiseuille equation, or the motion of a fluid through a porous media, 
described by the Darcy’s law (empirical law obtained by Darcy), its general solution 
remains one of the greatest challenges of mathematicians today. Therefore, the objective 
of this chapter is to present the resolution of the Navier-Stokes equation with the laws 
of wetting for different cases of interest in the manufacture of composite materials. 


Keywords: Wettability, Infiltration, Navier-Stokes 


1. Introduction 


Composite materials are an important and oftentimes used type of engineering materials with 
tuneable properties which depend on the materials combined. They are usually made of two 
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or more materials. The first phase, commonly referred to as “reinforcement” is constituted by 
a solid material, which can be employed in several morphologies such as fibers (short or 
continuous), particles, or as a co-continuous domain of a solid phase which occupies a fraction 
of the available volume (foams) [1-7]. The second phase, normally called “matrix,” can be 
gaseous, liquid, or solid, and is in contact with the reinforcement through a proper interface. 
Composite materials are very common in nature, where a combination of fibers and matrix is 
often found. A good example is wood, which is composed of flexible fibers of cellulose 
embedded in a rigid matrix of lignin. Artificially produced composites are usually classified, 
according to the nature of their matrix, into metal, ceramic, or polymer matrix composites. 
Each type of composite presents its advantages and disadvantages, and together they cover a 
great scope of applications [1, 8-12]. 


Composite materials manufacturing can be accomplished through several techniques. Each 
technique is mainly the result of the trade-off between the target properties of the final product, 
the starting materials properties, and the final cost. In most of the cases, the matrix is in a liquid 
state during processing. A common action in composite materials manufacturing is “to 
infiltrate” the reinforcement with the matrix. When a liquid is used, the impregnation of the 
remaining free volume is highly dependent on how the first phase is placed in the space and 
on the physical properties of the liquid itself. In the case of polymeric matrix composites (PMC) 
techniques involving a liquid matrix are the only way to process them, as polymers are in a 
liquid state before their cross-linking in thermosets, or can be melted if they are thermoplastics. 
Also, ceramic matrix composites are processed through liquid routes (e.g., sol-gel techniques 
for oxide composites, or polymer infiltration and pyrolysis (PIP) [13-14] for ceramic matrix 
composites (CMC) produced by pre-ceramic precursors. In the latter case, where several 
impregnations are performed, the permeability of the porous materials decreases with the 
increasing of the final density of the composite under production). A widely utilized technique 
which stands among ceramic and metallic matrix composites (MMC) manufacturing is the so- 
called reactive silicon infiltration [15]. Its peculiarity is that, during processing, silicon is molten 
and by reacting with extra carbon present in the primary phase it yields silicon carbide. This 
method was not originally conceived to produce composite materials, but for the synthesis of 
bulk SiC, still, it proved to be a suitable methodology for ceramic matrix composites manu- 
facture as well. 


Some industrial products obtained with this infiltration technology are pantographs, connect- 
ing rods, loudspeakers, piston rings, heat sinks, electro-brushes, engine cylinders, brake discs, 
etc. In most of these examples, primary phase (reinforcement) content is above 50% in volume. 
For reinforcement contents of such magnitude, the most commonly used technique for 
manufacturing composites is infiltration [1, 5-7, 16-18]. 


Actually, infiltration processes are not exclusive of materials manufacture, but can be observed 
in many everyday life activities, such as making coffee, where hot water is infiltrated into a 
porous bed of coffee particles. Perhaps one of the most important cases where infiltration 
phenomenon is observed is in the penetration of rainwater into the ground, which ends up 
filling aquifers that guarantee water supply to population in dry regions where the surface 
runoff is scarce. Because of the great importance of infiltration in geological sciences, it was in 
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this field where the deepest and most detailed studies of this phenomenon were performed. 
Many advances in this area of knowledge have been expanded and applied to other scientific 
fields such as materials science. The clearest example is Darcy’s law which is widely used in 
the design of composite materials. This empirical law was obtained by Henry Darcy whose 
original work published in 1856 is entitled: “Les Fontaines Publiques de la Ville de Dijon” (The 
Public Fountains of the City of Dijon). However, Darcy's law is not the only example that can 
be found in the literature. The relationship between water saturation and organic compounds 
in soils of different origins has resulted in the semiempirical equations of Brooks and Corey 
[19] and Van Genuchten [20], both widely used in composite materials. 


The aim of the present chapter is to introduce the reader to the basic concepts of wetting and 
infiltration, and to combine these concepts to solve the equations governing the infiltration 
process in several particular cases that often arise in the composite materials manufacturing 
industry. 


2. Experimental 


Although it will be discussed in detail in the upcoming sections of the chapter, it is important 
to clarify that infiltration, in most composite materials manufacturing processes, does not occur 
spontaneously (in the case of PMCs it does proceed spontaneously, but it is so slow that it is 
necessary to force it), so the liquid must be forced to infiltrate by applying an external pressure. 
From an engineering point of view, mainly two technologies have been developed for this 
purpose: the so-called gas-assisted pressure infiltration (GPI), where a gas is used to push the 
liquid; and another called squeeze-casting (SC) in which pressure is applied by means of a 
mechanical device. Each technique has its upsides and downsides, which will be discussed 
throughout this chapter. The main advantage of GPI over SC is that it allows a better control 
of physical parameters (pressure and temperature), and admits higher complexity in the 
manufactured parts. On the other hand, the main advantage of SC over GPI is that higher 
pressures can be achieved and the automatization of the process becomes easier. 


First the GPI technique will be analyzed, and for this purpose, the example of a laboratory 
equipment developed at the University of Alicante will be used [10, 11, 21-29]. A scheme of 
the equipment is shown in Figure 1. The equipment consists of a pressure chamber (maximum 
admissible pressure at 300 °C of 5 MPa) equipped with a resistance furnace, which is thermally 
isolated from the chamber by the use of conventional refractory material (temperature of 
chamber wall is always below 100°C). The furnace has a power of 1500 W and can reach 
temperatures of 1000°C. For measurement and control of furnace temperature, K type 
thermocouples, connected to a temperature controller, are used. The pressure control is 
entirely performed by a system of pneumatic and electronic valves connected to a pressure 
controller. The pressure vessel is made entirely of stainless steel and consists of a cylindrical 
body with an internal diameter of 30 cm, an outer diameter of 34.6 cm, a height of 45 cm, and 
a wall thickness of 3 cm. The closing of the camera is made with 12 M16 mounting studs. Both, 
the electrical feed (to thermocouple input and power cables) and the gas inlet and outlet for 
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pressurizing the chamber are positioned in the body. Additionally, the camera is connected to 
a primary vacuum pump, which allows the process to be performed in the absence of oxygen 
(low partial pressure). 


The manufacturing process of a composite material using GPI can be described as follows: a 
preform, which consists of a porous bed of compacted particles, is placed into a quartz (or BN- 
coated graphite) crucible and a metal ingot is placed on top of it (see Figure 1b). The crucible 
is then introduced inside the pressure chamber and vacuum is applied during at least 30 min., 
until a pressure of about 100 Pa is reached. The chamber is heated under vacuum at a very low 
rate (3°C/min) up to 250°C, to remove humidity and gas adsorbed on particles. Subsequently, 
further heating is applied up to a temperature 50°C higher than the melting point of the metal 
at a rate of approximately 5°C/min. The maximum temperature is then maintained for 30 min 
followed by argon injection to reach the desired infiltration pressure. Once the preform has 
been infiltrated (approximately 2 min. under pressure) the crucible is pushed towards the 
bottom of the chamber, which is cold, and directional solidification occurs, ensuring a high- 
quality composite (null porosity). Figure 2 shows a composite material obtained at the 
University of Alicante with this process, where it is possible to observe the absence of porosity 
and an intimate contact between reinforcement and matrix. 


Gas or vacuum inlet 
Safety valve 


Sample thermocouple Gas outlet 


Chamber thermocouple 
Sample crucible Graphite crucible 
Resistance furnace Metal 


Porous sample 


Detail of sample crucible 


Figure 1. Schematic representation of a gas-assisted pressure infiltration (GPI) system developed at the University of 
Alicante and arrangement of the materials (reinforcement and matrix) inside the crucible. 


At the University of Alicante, another equipment, in this case for squeeze-casting technique 
[30], was also developed. Figure 3 shows a schematic representation of the main body of this 
equipment and how it operates. The main feature of this apparatus is that it can work up to 
450°C and at a maximum pressure of 150 MPa. To resist these operating conditions, special 
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(W720 Bolher) steel was used to build the body of the equipment. Embedded resistances in the 
central body heat the system, and temperature is measured and controlled by K type thermo- 
couples connected to a temperature controller. 


10um 


Figure 2. Polarized light optical microscopy image of a composite material obtained by GPI. The composite material 
consists of needle coke (reinforcement) and copper (matrix). 
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Figure 3. Schematic representation of a squeeze-casting (SC) system developed at the University of Alicante. The oper- 
ating mode is indicated in the figure. 
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The manufacturing process of a composite material using SC can be described as follows: The 
system is first preheated to the operating temperature (between 350 and 450 °C depending on 
the alloy) to avoid early solidification of the metal matrix during infiltration. The preform is 
preheated (600 to 700°C) and the alloy is melted (150°C above the melting point) in two 
auxiliary furnaces. Thereafter, the preform is transferred to the main chamber, and the liquid 
metal is added. The desired pressure is applied as indicated in Figure 3. After infiltration, the 
system is cooled down and the manufactured material is removed as indicated in Figure 3. 


The main problem with this technique is that the processing is complex because the times and 
temperatures must be very accurate to keep the metal from solidifying before penetrating into 
the preform. This disadvantage is also an advantage, since the high velocity of the process 
generates a finer microstructure and thus improves the mechanical properties of the final 
material. 


Most importantly, it is possible to find in the market, both SC and GPI fully automated devices, 
but the basic operating principles are the same for the ones developed at University of Alicante 
with a scale factor between 2 and 10. 


3. Wetting basic laws 


One of the most relevant characteristics for the production of MMC through liquid processing 
routes is the wetting of the reinforcement by the liquid metals, which is quantified by meas- 
uring the contact angle between the reinforcement and matrix [31]. The contact angle between 
a flat substrate and a drop of liquid metal can be routinely measured by the method of the 
sessile drop. However, this measurement is not possible when the solid is in particulate form, 
in which case the wetting can be studied by means of liquid metal infiltration into compacted 
particulate samples [10, 11, 21-29]. 


Considering a solid, flat, non-deformable and chemically homogeneous surface in contact with 
a nonreactive liquid in the presence of a gaseous phase, if the liquid does not completely cover 
the solid, the liquid will form with the solid a contact angle 0, as depicted in Figure 4. 


Vapour 


Figure 4. Configuration of a sessile drop on a substrate within which it forms a contact angle 0. 
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The equilibrium value of O, which is used to define the wetting behavior of a fluid, can be 
written in terms of the balance of three energies present in the interfacial vapor-liquid-solid 
system, according to the well-known Young equation (1805): 


cos = 22 Ost (1) 


where g; is the interfacial energy between the phases i and j (s, l, and v represent the solid, 
liquid, and vapor phases, respectively). A contact angle of less than 90° identifies a liquid that 
wets the substrate, while larger angles identify liquids that do not wet it. An ideal wetting 
liquid will have an assigned value of O equal to 0°. 


The interfacial energy ow also called surface tension (y„), can be measured with sufficient 
accuracy by a large number of techniques, e.g., maximum bubble pressure [32, 33], pendant 
drop [34], or great crucible [35, 36]. However, the surface energy of the solid-vapor phase (o,,) 
and the solid-liquid (o4) phase can only be estimated. Consequently, the two variables that are 
unknown in the Young’s equation can be grouped into only one, called work of adhesion (W,), 
which expresses the strength of the solid-liquid bonds and was introduced by Dupré (1869): 


W, = Ow + Ov z Og (2) 


If now this equation is inserted in the Young’s equation, the so-called Young-Dupré equation 
is obtained: 


W, = 7p (1+ cos@) (3) 


All infiltration processes are driven by a differential pressure that allows the liquid matrix to 
penetrate into the porous reinforcement. This difference of pressure between the two fluids 
(liquid matrix and gas phase filling porosity) is defined as the capillary pressure (P.) and 


described by the famous Young-Laplace equation for a straight circular capillary: 


2-7, -cosO 
Pape (4) 


r 


Where y, is the liquid-vapour surface tension, 0 the contact angle, and r the capillary radius. 


Figure 5 schematizes different infiltration situations based on the wetting or non-wetting 
behavior of the system. 
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Wetting 0 << 90° Wetting 0 > 90° 


APc>0O == Spontaneous APc<0 = External Pressure 


Figure 5. Scheme of infiltration process for different wetting behaving systems. 


The energy shift that occurs when immersing the capillary into a liquid can be described in 
terms of work as: 


W, = Og = O (5) 


where W, is the work of immersion. For wetting systems (W;<0), infiltration proceeds spon- 
taneously without the necessity of applying an external pressure, as is the case of infiltration 
of silicon into carbon in the RBSC method. For non-wetting systems (W,>0), infiltration does 
not occur unless an external pressure, superior to the threshold pressure, is applied. The 
threshold pressure can be written in terms of work of immersion as: 


P, =S,-W, (6) 


0 i i 


where S; is the particle surface area per unit of volume of liquid matrix. Since 


Og) =.) Yp ' c0s0, The threshold pressure can be related to the contact angle as: 
P =5;*7,,°COSO (7) 


Since the measurement of the surface area of the particles becomes rather inaccurate when the 
particle diameter surpasses a certain value, and the surface roughness is very difficult to 
estimate, it is possible to assume a particular geometry for the porous media, e.g., cylindrical 
pores with a radius r as an approximation, in which case the surface area can be written as 


S;= eke and the threshold pressure can be written as: 
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iy op) 
P, =2-A-7,,cos@ rs (8) 


Where A is a geometrical factor introduced to describe deviation from the assumed geometry 
of the pores and the surface roughness of the material (it usually oscillates between 2 and 4 [10, 
21-23, 25, 29], y,, is the surface tension, ¢ the porosity of the sample, r the mean pore radius, 
and 0 is the contact angle at the triple line. It is important to emphasize that the contact angle 
in eq. 5 is the equilibrium contact angle that can be measured with high accuracy in sessile 
drop experiments. 


4. The simplest case: Obtaining Darcy’s law 


As previously mentioned, Darcy’s law predicts how a fluid moves through a porous system. 
This empirical law is based on the study of water movement in sand beds, and represents an 
ideal model in the manufacture of a composite material. Darcy’s law can be derived from the 
resolution of the Navier-Stokes equation, considering the ideality of the system. In this section, 
the basis for describing the process of manufacturing of composite materials will be provided, 
and the most important factors that affect the process will be identified. 


As aforementioned, the flow of viscous liquids is governed by the Navier-Stokes equation. 
This differential equation does not have a general solution. However, in some cases, it is 
possible to solve, e.g., for an incompressible fluid, in laminar flow (low Reynolds number), 
assuming that the time derivative of the fluid velocity is much smaller than its spatial deriv- 
atives, and that the effect of gravity is negligible. In this case, the equation can be reduced to: 


VP =7-V°v (9) 


where v is the speed of the fluid, n is the dynamic viscosity, and P the pressure. Assuming 
unidirectional flow (which is not so realistic in a standard porous system), it is possible to 
reduce the expression to the Darcy’s law. This equation can be rewritten in the most commonly 
used form to analyze the flow of incompressible fluids through porous media: 


= (10) 


where z is the flow direction, v, the average superficial velocity of the fluid, dP/dz the pressure 


gradient along the infiltration front, and K the permeability of the porous media. The super- 
ficial velocity v, can be converted into infiltration velocity (dz/dt) in the porous medium by 


means of the porosity ¢: 
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V =o (11) 


Combining both equations and integrating, the Darcy’s law, which shows the relationship 
between infiltration distance (h), time (t) and the pressure drop in the infiltrating liquid (AP=P- 
Po) is obtained. 


h? = 2Kt Ap (12) 


1 


It must be pointed out that the permeability of the porous media is proportional to the square 
of the average pore radius: 


K= ar? (13) 


Darcy’s law has been widely used in basic studies of manufacture of metal matrix composites. 
Despite its simplicity, it has proven to be very useful in metallurgical engineering. As an 
example, the results of infiltration of liquid aluminium into a porous bed of particles, obtained 
by Narciso et al., are shown in Figure 6 [10, 25]. Figure 6a provides the relationship between 
infiltration distance and applied pressure for liquid aluminium into a porous bed of alumina 
particles, while Figure 6b shows the results for a liquid infiltration of an aluminium alloy into 
preforms of graphite particles. As predicted by Darcy’s law, the square of infiltration distance 
is proportional to the pressure drop (h? vs AP) in these systems. 


h (mm) 


1000 1040 1080 1120 1160 0 20 40 60 80 100 
P (kPa) AP(kPa) 


(a) (b) 


Figure 6. Evolution of the infiltrated height as a function of pressure for liquid aluminium on a preform of alumina 
particles (Alumina particles AA10, Vp = 59 %, Temperature = 700 °C) (a). Evolution of the square of the infiltrated 
height for a liquid aluminium alloy as a function of pressure drop on a preform of carbon particles at temperature of 
700 °C (see the appendix for the characteristics of the graphite particles) (b). 
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As discussed in the preceding paragraphs, to solve the Navier-Stokes equation a full saturation 
and a laminar flow (Re < 2100) have been assumed. Besides, the effect of gravity has been 
considered to be negligible. Figure 6b shows a good correlation of the experimental data with 
Darcy’s law, and thus, the assumed conditions should be valid. Reynolds number (Re) can be 
obtained from the following expression: 


_ pu:D 
n 


Re (14) 


Where p is the density of the liquid metal; D is the diameter of the channel, which in this case 
is the interparticle distance that can be assumed to be about 1/3 of the average size of the 
particles; v is the fluid velocity, which is comprised between 1 and 5 mm/s (in our equipments); 
and n the dynamic fluid viscosity. If the Reynolds number is calculated for the experiments 
shown in Figure 6, the values obtained are between 0.001 and 0.7, which means that the laminar 
flow assumption is indeed correct. The force exerted by gravity would be the metallostatic 
pressure of the melt, which is given by: 


Papen (15) 


And considering that the liquid column is always less than 10 cm, in the case of aluminium 
for example, the pressure would be about 2 kPa, which is almost 3 orders of magnitude lower 
than the capillary forces, and therefore negligible. 


In the following sections different real situations that can occur are considered and will be 
analyzed separately. The combination of all these particular cases simultaneously, although 
more realistic, is much more complex and goes beyond the scope of this chapter. 


5. The effect of saturation 


On systems where the liquid does not wet the solid, which is most common in MMCs, it must 
be taken into account that the permeability of the porous system is not constant [37-39] and 
that it depends on the saturation and this latter, and at the same time, on the applied pressure 
[21, 40-42]. 


The study of unsaturated flow is well known in geology and in 1931 it was described mathe- 
matically by Richards [34], whose partial differential equations are difficult to solve and 
usually solved numerically for each particular case. Numerous researchers have worked on 
this particular hydrodynamics and mathematical problem [44,45], and the books by Prof. Bear 
collect these study cases [46-48], which are essential reading for researchers working in this 
field. It is also noteworthy that Prof. Bear highlights in his books the goodness of Darcy’s law, 
and that its application range is wider (more useful) than it was expected at first (i.e., unsatu- 
rated flow, high Reynolds number, etc.). 
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These studies were not limited to the case study of hydrology, but have crossed over to the 
field of metal matrix composites (MMC), performed especially by the group of Andreas 
Mortensen [40-42], but as mentioned above without analytical solution of the differential 
equations describing the phenomenon. The highlight to be found in the recent literature are 
two articles that show that close to the threshold pressure, infiltration is actually governed by 
percolation phenomena (instead of slug flow) [49,50]. 


Another aspect to consider is that pressure also depends on the position along the flow path, 
and is usually assumed that the pressure drop between the front (P,) and head (applied 
pressure) is linear, which is reasonably sane. Nevertheless, in order for the differential 
equations to have an analytic solution, it will be assumed in this chapter that practically all the 
preform is at the same pressure, which in principle is something reasonable as working 
pressures oscillate between 1 and 5% of the threshold pressure (i.e., head=1.05P,; front=P,). 


In the literature, several models showing the relationship between pressure and saturation 
have been described, but the most famous are the Brooks-Corey (BC) [19] and van Genuchten 
(VG) [20] models. In this section, the Navier-Stokes equation has been solved considering that 
the permeability is not constant, relating it to the saturation through the BC and VG models. 
Different authors have used these models to describe the infiltration process in MMCs [51]. 
However, previous works have shown that for low pressures (close to the threshold pressure) 
these models, especially the BC, do not describe the infiltration process correctly [52, 53]. 


Starting from reduced Navier-Stokes equation: 
dz? = ——<dP (16) 


Where z is the infiltration distance, K the permeability, t the time of infiltration, 7 the dynamic 
viscosity of the fluid, Ø the porosity of the porous body, and P the applied pressure. Perme- 
ability can be decomposed into two different terms, intrinsic permeability (K;), which depends 
exclusively on the nature of the solid phase; and relative permeability (K,), which depends on 
the degree of saturation (S) of the system: 


K=K;,K, (17) 


1 r 


Relative permeability, depends on the degree of saturation according to: 


K = s" (18) 


r 


Where n is a parameter that ranges from 1 to 3. In the case of composites fabricated by 
infiltration, it has been found that 1 is the most accurate value for n [51]. There are several 
models that relate saturation to the applied pressure (P) and the threshold pressure (Po) of the 
system. One of the most used models of this kind is the Brooks-Corey equation: 
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S= 1-[2) (19) 


Where A is a parameter that depends on the solid material. Replacing K in eq. 16 with eqs. 
17, 18, and 19 yields a modified reduced Navier-Stokes equation: 


p y i 
«1-(3] | i (20) 
=— ` ^d gp 
URA 


Solving this equation provides a model to estimate depth of infiltration as a function of applied 
pressure. In this form, there is no general analytical solution to the equation, so, to obtain an 
analytical model, a particular case must be considered. As discussed above, experimental data 
fit better when n=1, so to solve the equation, this particular case will be considered: 


a 
1 K,+t P 
Zh =i |} 1-} 2 -dP zj 
2 na j (2) | en 
a 
{2 
aK P 
yea ON) py 7 a (22) 
na A-1 


Where A is the integration constant, and its value can be determined using initial conditions, 
i.e., for P= P, and h=0 


wa) A (23) 
0=P +—— +A 
P 
A= (24) 
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Then the final expression for h(P) is: 
{2 
2-K;t P 1 
k ==] P+— -P — +1 25 
h-rd) 05) 


It must be noted that if A acquires a very high value, the model tends to Darcy’s law: 


A >æ 
h? = FApo -P,(0+1)) 
m 
2-K;+t 
h? ee (P= 


If, instead of using the Brooks-Corey equation, the Van Genuchten model is applied: 


Şell- [1 +(a(P- py] (26) 


ata (27) 


where a is a scale parameter inversely proportional to the mean pore diameter, and 11 (called 
n1 to differentiate from n in eq. 18) and m are shape parameters. The reduced Navier-Stokes 
equation results in: 


ae -[1+0-(P-P Te] : (28) 
Pp 


URA 


dz? = 


As in the case of the Brooks-Corey model, a general analytical solution is not available, so the 
case of n = 1 will be considered again: 


el oes O 


URA 


dz? = 
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zi =] 1-fi+((p -p T" fap (30) 
Ly = ~ frap- f [i+ (aP -pe ap (81) 
zh = oP - [[1+(a-(P- py] dP (32) 


The second integral does not have an analytical solution either, so in order to obtain a solution 
as general as possible, the value of n1=2 is introduced [42, 51]: 


Ket -1/ 
2 1 *| == . — 
D [P f[t+(@(P-P,) | aP | (33) 
Lpa Atip f l dP (34) 
2 nD J1+(a(P-P,)P 
Applying the substitution u=a -(P-P ), du=a-dP results in: 
n =A pL j- ay (35) 
URA 1+4? @ 
2K: 
h? = SKE P= 1 1 du (36) 
12 a" /1+u? 
2.K arcsinh (u 
pa | (o) c (37) 
RA a 


Replacing back u=a-(P-P) : 
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K- inh(a-(P — P. 
Kt p arcsi (a-( D) 
RZ a 


h? 


€ (38) 


To estimate the value of C, initial conditions must be applied P =P, and h =0. 


2 2:-K;t ) arcsinh (0 ) 


0 Pr F T T A Ned (39) 
C=-P, (40) 


Then, the general solution for n1 = 2 is: 


pakti p-b) arcsinh(a-(P — P,) ) 
B ea a ae 


a a (41) 


In order to analyze the accuracy of the results, the evolution of the square of the infiltrated 
height as a function of pressure drop is plotted in Figures 7a and 7b. As shown, experimental 
results are perfectly framed between the line defined by Darcy’s law (assuming that the liquid 
permeability is invariant and therefore permeability measured with polyethylene glycol is 
used) and the curve modified with the Van Genuchten equation. 
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Figure 7. Evolution of the square of infiltrated height for liquid aluminium-silicon alloy as a function of pressure drop 
on a preform of graphite particles P1 at 700°C) (a). Evolution of the square of the infiltrated height for liquid alumini- 
um-silicon alloy as a function of pressure drop on a preform of carbon particles P2 at 700°C) (b) Characteristics of the 
graphite particles are listed in the appendix. The coefficients for the BC model have been extracted from [53]. 
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The results shown here clearly demonstrate that it is not necessary to solve the Richards 
equation (modified Darcy’s law to take into account saturation), which has only numerical 
solutions to predict infiltration in composite materials unlike what was suggested by Morten- 
sen et al. [54]. It is noteworthy that many authors in soil mechanics indicate that Darcy’s law 
is more universal than it was originally supposed. 


6. The effect of variable radius 


Although the variation of the radius is not a very common phenomenon in the manufacture 
of composite materials, it is a process that can occur. The variation of the radius is a funda- 
mental issue for example in the manufacture of SiC by the RBSC process or reactive infiltration 
and SiC derived CMCs. This is a rather complex situation since the variation of the radius 
affects the permeability of the porous media. 


To solve the Navier-Stokes equation it has been assumed that the reduction of the radius, 
caused by a chemical reaction, is controlled by diffusion. The equations presented so far are 
valid for infiltration in nonreactive systems. However, they are not valid for reactive infiltra- 
tion because the pore radius, in this kind of systems, is not constant but decreases with the 
generation of new solid. As discussed above, the threshold pressure and the permeability 
depend on the radius, so they are not invariant over time. 


So for reactive systems, it is necessary to re-solve the reduced Navier-Stokes equation consid- 
ering the variation of the radius: 


pe = (42) 
After integration, the following expression is obtained: 
K 
—h?=- ——dP-dt 43 
5 ah (43) 


Usually in reactive systems, infiltration takes place spontaneously (without applying an 
external pressure), so AP = P}. Differentiating the expression for P} results in: 


T= 
dP =-2-A-y,, coo Dy (44) 
T 


Combining this expression with eq. 6 yields: 
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1 ar? (1-4) 
h’ = 2-A-y,,-cos0- -dr-dt 45 
2 Qn | Vio Or ( ) 
Simplifying and rearranging: 
Ta 
k = MNA. 6), dr-dt (46) 
12 


In principle, porosity of the system is not constant, but for a first approximation it can be 
considered so. To proceed with the solution of the equation, a relationship between radius and 
time is necessary. Assuming that the controlling step of the reaction is the diffusion of C 
through the SiC layer, then the model of decreasing nucleus developed in the decade of the 
1950s can be used [55, 56]: 


2 3 
M.-0. r? 
[= si Psich t-32) (4) (47) 
6M ic. Deysic Psi Ti r, 
M-00.. r° 2 
P a Bae CAKES dr (48) 
6-Msic.Pojsic'Psi n n 


Where pṣ;c is the density of SiC, rọ the initial pore radius, D¢;sic the diffusion coefficient of C 
atoms through SiC, and ps; the density of Si atoms in kg/m*, Msic and M; the molar weight of 
SiC and Si, respectively. 


Introducing these relations into eq. 46 and defining: 


(1 ‘| 9) M;i Psic 1 
ne 6-De sic’ Mic Psi 


G=4-/A-a7,,:cos0- (49) 


Results in the following expression: 


h= G ka + of 5 ra (50) 
r r 


Solving the first integral yields: 
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h? = cj a-ofž] +22) dr (51) 


0 0 


Where A is an integration constant. Solving the second integral results in: 


>) of Pid (5) X) H5) A J (52) 


where B is an integration constant. In order to obtain the final solution, boundary conditions 
must be used to find the values of A and B: 


a r=r, ~h=0 -t=0 


b. dh 


r=0 -h=hyay ~t=t— G=0 


Ms; ` Psic’ p 
Where t= 
6- M sic. 2ysic * Csi 


By applying condition a), the next expression is obtained: 


ry?) 1 
j= ofan- a (53) 


Applying condition b) leads to: 
hmax = GB (54) 


and 


"T Gan a a S 
d 3 4 (55) 
Í 2 cf 2 8 $ 5] Fe J 
r 2r 


Solving eq. 55 leads to: 
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Replacing the value of A in eq. 53 results in: 


Replacing the value of B in eq. 54 provides the expression for hyay: 


(1 — 9) M;;Psic A 
n2? 3-Do sic Msic'Psi 


hax = | Aay ,,,COSO (56) 


And after replacing the values of A and B, and rearranging the final expression h(r) is obtained: 


1- we 3 2 
h? = E n Psic Ti r H r 2. r 4 1% (57) 
IRA 3°De sic’ Psi f % r 
42 3 2 
h2 = IMAX y. r 2. r i To (58) 
% % % f 


4 3 
A gl E EE EA 
AJA ‘ 


It must be taken into consideration that this equation is valid only as long as the modification 
of the porosity is negligible. For evaluating h(t), eq. 59 and eq. 47 must be combined. 


Figure 8 shows the calculated infiltration distance with time for a carbon preform with an initial 
pore radius of 30 micron and 50% open porosity, infiltrated with pure Si, using Darcy’s law 
and the variable radius model for cylindrical and spherical pores. The variable radius models 
predict a slower rate of infiltration, which is closer to the experimental data obtained by several 
authors [12, 57-61]. 


Another approach is to consider that the pore closure is limited by the chemical reaction that 
originates it. In this case, the equations of the decreasing nucleus model to be used are: 


= M3; Psicho A 1 = z (60) 
k Psi Msi 


iC 
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Figure 8. Comparison of the square distance of infiltration (h?) vs. time, calculated using Darcy’s law and variable radi- 
us models for a carbon with an initial pore radius of 30 micron, and 50% porosity, infiltrated with pure Si. For the 
variable radius model, the pore closure is assumed to be controlled by diffusion, and two different pore geometries are 
considered. Infiltration is simulated until the pore radius decreases by 30%. 


Psic 
dt = — —>=—-dr 61 
k, ‘Psi ( ) 


Introducing this equation into the reduced Navier-Stokes equation and solving the integrals 


results in: 


1- sie 
h’ = S a = FEE : va fe + A)-dr (62) 
7 3s Pg Vi sic 


1-¢) M,-p,; : 
k = ~ 07 co 9) Matas | + Art J (63) 


The integration constants A and B are calculated using contour conditions: 
ce r=r, ~h=0 >t=0 
d. r=0 >h=hmax 7 t=T>dh /dr=0 


a. For condition a): 


I= M.j. De. a 
0= E 9) MePa {n +Ar,+B (64) 
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For condition b): 


= 
eas o), Mai Por (0+4) 
dh Se n2 ke Ps Msic (65) 
dr J= . 
2- PE coed: 6), Ms'Pac (3 +A0+ B) 
WO” kPa Msc \2 
A=0 (66) 


and replacing in condition a): 


1-$¢) Mp; : 
0= P ca) 9) _MaPac Psic |4 O-r, + B (67) 
nZ Kp Msc \ 2 
2 
pss a (68) 
and 
1- ape. 
hax = fn cose! a Msi Psic y2 (69) 
7” ke ps Mic 
Replacing the values of A, B and hyax: 
1+ M d f 2 y2 
h? = tray eos I) Mar Pic mE (70) 
i n2 kePsi Msc 2 


Figure 9 shows the square of the infiltration distance as a function of time for a carbon preform 
with an initial pore radius of 30 micron and 50% open porosity infiltrated with pure Si, using 
Darcy’s law and the variable radius model for a chemically controlled pore closure. As seen, 
infiltration rate decreases with time, just as in the previous case (diffusion controlled pore 
closure). However, the experimental data [12, 42, 59-66] show that the infiltration rate is lower 
than the one predicted by these models; furthermore, infiltrated height is proportional to t 
instead of to t". Experiments indicate that infiltration is controlled by the reaction at the triple 
line and not by the viscous phenomena described here. However, both interpretations are 
compatible. It is probable that in the first stage of the process, infiltration is controlled by the 
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reaction at the triple line, and after the SiC interface has been formed, the viscous flow becomes 
the controlling process. A completely different approach to the process of infiltration is 
performed by Sergi et al. [67, 68], in which they use Lattice-Bolztman models to describe 
infiltration in a single capillary and predict that infiltration distance is proportional to time 
when considering pore narrowing of the capillary due to the reaction. However, their predic- 
tions have not been validated experimentally up to date. It is possible to conclude that a deeper 
understanding of the reactive infiltration phenomenon should be gained in order to discern 
what is happening in the overall process. 


0 2 50 75 100 125 150 175 200 


t (s) 


Figure 9. Comparison of the square distance of infiltration (h?) as a function of time, calculated using Darcy’s law and 
the variable radius models for a carbon with an initial pore radius of 30 micron, and 50% porosity, infiltrated with pure 
Si. For the variable radius model, pore closure is assumed to be controlled by the chemical reaction. Infiltration is si- 
mulated until the pore radius decreases by 30%. 


7. Effect of variable viscosity: Curing polymeric resins 


The properties of PMCs depend strongly on the type of matrix used. Among the possible 
binders, thermoset resins present a singularity: their viscosity is time dependent [69, 70]. These 
kinds of polymers present an initial low viscosity that grows exponentially with time. During 
this process, the material undergoes a series of chemical reactions that cause the cross-linking 
of polymeric chains yielding a 3D structure with excellent thermal stability and solvent 
resistance. The process of cross-linking that takes place in these resins is known as curing. In 
polymer matrix composites in which the matrix cures with time, thus hardening and losing 
fluidity, viscosity becomes a limiting factor in infiltration. 


The evolution of viscosity is usually expressed with the following equation: 


= peT va 
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Where ngis the initial viscosity and k, is a strongly temperature-dependent rheological kinetic 


constant [70]. Starting from the reduced Navier-Stokes equation: 
Q— =-—— (72) 
After integration, the following expression is obtained: 
1 K 
—W=- —dP-dt 73 
>. Ba (73) 


Introducing eq. 71 and solving the integral yields a model for predicting infiltration distance 
as a function of time for curing polymeric resins: 


ix; K 
7 h one aP-dt (74) 
0 
1 K dP 
A ee (75) 
0 
W=- Sa feat (76) 
0 
kt 
pao Sa nee (77) 
ON n 


To calculate the values of integration constants A and B, initial conditions are applied: 
a. t=0;h =0 


b. 428 2 
t=;  =0 


For condition a): 


(78) 
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for condition b): 


a _2KAP it, a 
7 
af pä (e2) 
A=0 


Then the final equation results in: 


pa 2K APs KAP 


= ——e ” — 
D-k, Dnk, a 
2-K-AP He 
h? s . 1 e k,, t 
Fae ( ) (81) 


It must be highlighted that the negative sign is only indicative of the direction of flow. Figure 
10 shows the evolution of the square of infiltrated height for a preform infiltrated with curing 
resins with an initial viscosity of 100 Pa s at different temperatures. The drastic effect of 
temperature on infiltration distance becomes evident, since increasing temperature to 50°C 
reduces infiltrated distance by one order of magnitude. 
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Figure 10. Evolution of the square of infiltration height with time for curing resins into porous preforms at different 
temperatures. 
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8. Concluding remarks 


The present chapter has introduced the reader to the processes involved in the manufacture 
of different types of composite materials (PMCs, MMCs, CMCs) via infiltration. Over 90% of 
the world production of composite materials is carried out by infiltration, especially in the case 
of PMCs where it is almost the only available processing technique, while in case of metal 
matrix, production via infiltration represents only 70% [7]. In order to accurately control the 
properties of the final material and to be able to optimize the production, it is key to identify 
and understand the different processes involved in the manufacture of these materials. For 
this purpose, in this chapter, different models have been elucidated, starting from the differ- 
ential equation of Navier-Stokes and taking into account the basic laws of wetting, by applying 
different real conditions often found in industry. For PMCs, the main discussed issue was the 
change in viscosity that thermoset matrixes experiment during infiltration. By introducing this 
phenomenon in the infiltration equations, it becomes clear that it can be a limiting factor during 
production since it shows a great deviation from Darcy’s law. In the case of MMCs, the main 
phenomenon analyzed was the effect of saturation on infiltration. It became clear that variable 
saturation not only affects in a huge manner infiltration kinetics, but also changes the domi- 
nating phenomena. The last publications in this area have identified percolation as probably 
the leading phenomenon on infiltration in unsaturated systems at very low overpressure. 
Finally, for CMCs, the role of pore narrowing, an issue that has been widely discussed in the 
literature over the past 10 years was revised. This phenomenon is especially relevant in 
processing of ceramics and ceramic matrix composites via reactive techniques like the RBSC 
method used for SiC production. Even though the research on this area is vast, there is still not 
a full understanding of the effect of pore reduction on infiltration. Most models proposed still 
yield higher infiltration rates than the ones observed experimentally, because they do not 
consider the effect of chemical reactivity. 


Appendix 


Table 1 shows the main characteristic of the carbon particles used in this work. 


Particle D(90) D(50) D(10) D(4,3) D@G,2) Span V, 
P1 25.8 13.4 7.0 15.1 11.6 1.39 0.511 
P2 41.0 25.8 15.5 27.2 23.0 0.99 0.515 
P3 95.1 62.0 36.8 64.0 53.0 0.94 0.518 
P4 173.6 117.6 56.6 124.1 110.4 0.99 0.482 


Table 1. Main characteristics of the carbon particles used. Span is defined as D(90) - D(10)/D(50), where D(x) 
represents the diameter below which x% of particles are encountered. V, is the particle volume fraction. 


Table 2 shows the value of the parameters required to assess h(t) for the several cases exposed 
in this work. All parameters are in SI units. 
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Symbol Magnitude Value SI units Ref. 
A Geometrical factor 2-4 - [10, 21-23,25, 29] 
a Proportionality constant 4-104 [61] 

Yr Liquid-vapor surface tension of Si (at melting point) 0.750 N/m [34-36] 
0 Contact angle of Si/C (at Si melting point) 40 2 [31] 
n Dynamic viscosity of molten Si (at melting point) 0.605-10° Pa's [34] 

Qsic Density of SiC 3210 kg/m? [9, 12] 

Mic Molar weight of SiC 40-10% kg/mol [9, 12] 
Ms Molar weight of Si 28-10% kg/mol [9, 12] 
csie Diffusion coefficient of C through SiC 2-10% m?/s [62] 

Ost Density of Silicon (at melting point) 2570 kg/m? [65] 

k kinetic constant in the SiC formation reaction 4-108 m/s [71, 72] 


Table 2. Value of the main parameters used in this work in SI units. 
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